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GAS  MOTION  IN  A  LOCAL  SUPERSONIC  REGION  AND  CONDITIONS 
OF  POTENTIAL-FLOW  BREAKDOWN* 
By  A.  A.  Nikolskii  and  G.  I.  Taganov 

For  a  certain  Mach  number  of  the  oncoming  flow,  the  local 
velocity  first  reaches  the  value  of  the  local  velocity  of  sound 
(M  =  1)  at  some  point  on  the  surface  of  the  body  located  within 
the  flow.  This  Mach  number  is  designated  the  critical  Mach  number 
M(,j..  By  increasing  the  flow  velocity,  a  supersonic  local  region 

is  formed  bounded  by  the  body  contour  and  the  line  of  transition 
from  subsonic  to  supersonic  velocity.  As  is  shown  by  observations 
with  the  Toepler  apparatus,  at  a  certain  flow  Mach  number  M  >  M^j., 

a  shock  wave  is  formed  near  the  body  that  closes  the  local  super- 
sonic region  from  behind.  The  formation  of  the  shock  wave  is 
associated  with  the  appearance  of  an  additional  resistance  defined 
as  the  wave  drag. 

In  this  paper,  certain  features  are  described  of  the  flow  in 
the  local  supersonic  region,  which  is  bounded  by  the  contour  of  the 
body  and  the  transition  line,  and  conditions  are  sought  for  which 
the  potential  flow  with  the  local  supersonic  region  becomes  impos- 
sible and  a  shock  wave  occurs. 

In  the  first  part  of  the  paper,  the  general  properties  of  the 
potential  flow  in  the  local  supersonic  region,  boionded  by  the  con- 
tour of  the  profile  and  the  transition  line,  are  established.  It 
is  found  that  at  the  transition  line,  if  it  is  not  a  line  of  dis- 
continuity, the  law  of  monotonic  variation  of  the  angle  of  inclina- 
tion of  the  velocity  vector  holds  (monotonic  law).  An  approxima- 
tion is  given  for  the  change  in  velocity  at  the  contour  of  the  body. 
The  flow  about  a  contour  having  a  straight  part  is  studied. 

In  the  second  part  of  the  paper,  an  approximation  is  given  of 
the  magnitudes  of  the  accelerations  at  the  interior  points  of  the 
supersonic  region.  With  the  aid  of  these  approximations,  it  is 
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shown  that  for  profiles  convex  to  the  flow  the  breakdown  of  the 
potential  flow,  associated  with  an  increase  of  the  t-Iach  munber  of 
the  oncoming  flow,  cannot  be  due  to  the  formation  of  an  envelope 
of  the  characteristics  within  the  supersonic  region. 

On  the  basis  of  the  monotonic  law,  the  transitional  Mach 
number  M  is  found,  beyond  which  the  potential  flow  with  local 
supersonic  region  becomes  impossible. 


I  -  PROPERTIES  OF  POTENTIAL  FLOW  WITH  LOCAL  SUPERSONIC  REGION 

1.  Properties  of  Supersonic  Flow  Bo\inded  by  Solid  Wall 

and  Transition  Line 

The  flow  in  the  supersonic  region  about  a  wing  profile  of 
usual  form  with  a  flow  about  it  of  a  Mach  number  higher  than 
critical  is  considered.   In  accordance  with  these  conditions,  two 
types  of  such  flow  may  occur: 

(a)  Transition  from  supersonic  flow  to  subsonic  flow.  -  In 
this  case,  the  supersonic  region  is  bounded  by  part  of  the  solid 
contour  AB  and  the  transition  line   X  =  1,  which  is  indicated 
by  a  dotted  line,  as  shown  in  figure  1.  The  magnitude   X  =  w/a^, 
where  w  is  the  absolute  magnitude  of  the  velocity  and  a*  is 
the  critical  velocity  (the  velocity  at  which  the  flow  velocity  is 
equal  to  that  of  sound); 

(b)  Transition  from  supersonic  flow  to  subsonic  flow  with 

aid  of  shock  wave.  -  In  this  case,  the  supersonic  region  is  bounded 
by  the  solid  contour  AB,  the  transition  line  X  =  1,  and  the 
shock  wave  CB  (fig.  2). 

Through  each  point  of  the  supersonic  region,  two  character- 
istics of  different  families  pass.  The  minimum  angle  between  the 
direction  of  the  velocity  at  a  given  point  and  each  of  the  char- 
acteristics is  equal  to  the  Mach  angle  a. 

If  the  velocity  vector  is  rotated  by  the  angle  a  in  a 
counterclockwise  direction,  the  direction  of  the  vector  will  coin- 
cide with  one  of  the  characteristics  that  shall  be  denoted  as  the 
characteristic  of  the  first  family  in  contradistinction  to  the 
other  characteristic,  which  shall  be  denoted  as  the  characteristic 
of  the  second  family. 
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In  case  (a),  each  of  the  characteristics  drawn  from  an  arbi- 
trary point  P  of  the  contour  has  a  point  in  common  with  the 
transition  line  X  =  1  (fig.  1). 

In  case  (b),  all  the  characteristics  of  the  second  family  end 
on  the  transition  line,  but  not  all  characteristics  of  the  first 
family  possess  this  property. 

By  displacing  the  point  P  on  the  contour  in  the  direction 
toward  the  base  of  the  shock  wave  B,  it  is  fo\ind  that  originating 
from  a  certain  point  D  of  the  contour,  the  characteristics  of 
the  first  family  no  longer  end  on  the  transition  line  but  intersect 
the  shock  wave  (fig.  2).  The  essential  difference  of  the  flow  with 
local  supersonic  region  from  other  mixed  flows  (for  example,  flow 
from  a  Laval  nozzle)  lies  in  the  fact  that  from  each  point  of  the 
contour  there  originates  at  least  one  characteristic  that  ends  on 
the  trains  it  ion  line. 

The  region  ADCA,  in  which  both  characteristics  drawn  from  the 
points  of  the  contour  end  on  the  transition  line,  is  denoted  as 
region  I;  the  region  DCBD,  in  which  only  the  characteristics  of 
the  second  family  fall  on  the  transition  line,  is  denoted  as 
region  II. 

In  case  (a),  the  entire  supersonic  region  coincides  with 
region  I,  hence,  all  theorems  derived  in  thxs  paper  and  the 
approximations  for  the  flow  in  region  I,  obtained  from  the  single 
assumption  of  the  ending  of  the  characteristics  of  both  families 
on  the  transition  line,  will  also  hold  for  case  (a). 

Hereinafter,  there  will  often  be  considered,  together  with  a 
certain  point  P  of  the  contour,  simultaneous  points  on  the 
line  X  =  1,  which  are  the  ends  of  the  characteristics  originat- 
ing from  point  P.  The  end  of  the  characteristic  of  the  first 
family  will  always  be  designated  by  the  same  letter  as  the  point 
of  the  contour  but  with  subscript  1  and  by  the  sign  *.  Similarly, 
the  end  of  the  second  characteristic  will  be  designated  by  sub- 
script 2  with  sign  *,  as  shown  in  figures  1  and  2. 

In  the  following  derivations,  the  fundamental  magnitudes  are 
the  inclinations  of  the  tangents  at  an  arbitrary  point  P  of  the 
contour  to  the  axis  of  abscissas,  denoted  always  by  0^;  the 

inclinations  of  the  velocity  vector  on  the  transition  line  at 
points  P,   and  Pg^  are  denoted  by  9-^^     and  Qg*^  respectively. 

All  three  magnitudes  9-^,     0,^,  and  Sg*  ^^®  functions  of  the 
arc  length  of  the  contour. 
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An  arbitrary  point  P  on  the  contour  in  the  supersonic  region 
is  considered.  The  characteristic  of  the  first  family  PPg*  is 
represented  in  the  plane  of  the  hodograph  of  the  velocity  by  an  arc 
of  an  epicycloid  of  the  second  family  P'Pg*'  (fig.  3).  Let 

6  =  02^^  -  ^v  ^®  ^^^   polar  angle  between  the  points  Pg**  ^^^ 
P]^^'.  For  the  magnitude  X  at  the  point  P  of  the  profile,  the 
following  equation  is  obtained: 

X  =  f{e^  -  e^)  (1-1) 

where  X  =  f (6)  is  the  equation  of  the  epicycloid  so  set  up  that 
the  equation  f(0)  =1  holds.  This  equation,  as  is  known,  has  the 
form  (reference  2) 


e  =  _  ± 


i  {\IW  ^"''^  ^^"^  & "  ^^"^^  ^^J 


arc  sin 


"'^^'  r^-^)<JB-^)l    f^-^) 


In  determining     Sg*     from  the  given     X    sind     9-^ 

^2*  =  ^k  +  <^M  (1.3) 

where  cp  is  the  function  reciprocal  to  f .  The  graphs  of  the 
functions  f(9)  and  cp(X)  are  given  in  figure  4. 

By  considering  the  characteristic  of  the  first  family  starting 
from  point  P  (fig.  3), 

X  =  f  (ej^  -  e^^)       01^  =  0j^  -  cp(X)         (1.4) 

where  f  and  cp  are  the  ssime  functions  as  in  equations  (l.l)  and 
(1.3). 

For  the  points  of  region  I  where  the  characteristics  of  both 
families  end  on  the  transition  line,  the  following  equation  is 
obtained  from  equations  (1.3)  and  (1.4): 
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Thus 


TBEOEEM  1.  -  The  inclination  of  the  contour  at  any  point  F 
of  zone  I  is  the  arithmetic  mean  of  the  inclinations  of  the  veloc- 
ity vectors  at  the  points  of  the  transition  line^  which  are  the 
ends  of  the  characteristics  starting  out  from  point  P. 

Equation  (1.3)  and  the  second  equation  (1.4)  permit,  for  a 
certain  distribution  of  the  velocities  on  the  given  contour,  deter- 
mination of  the  inclinations  of  the  velocity  vectors  at  the  ends 
of  the  characteristics  lying  on  the  transition  line. 


2.  Monotonia  Law  of  Change  in  Angle  of  Inclination 

of  Velocity  Vector  on  Transition  Line 

If  in  the  equations  of  an  adiabatic  gas,  the  nondimens ional 
velocity  X  and  the  angle  of  inclination  9     of  the  velocity 
vector  to  the  x  axis  are  taken  as  the  xxnknown  fvinctions,  the 
following  equations  are  obtained: 

(1  -  m2)(co6  e|^  +  sin  e^  -    xLn  6^  -   cos  9^)  =  0 

sin  9^  -  COB  ^  +   X  cos  0^  +  \  sin  e^^  =  0       (1.6) 
ox       oy        dx        oy 

The  character  of  the  change  of  the  magnitude  '  9     along  the 
line  X  =  Xj^  =  constant.  On  one  side  of  this  line,  let  A  <  X]^ 

and  on  the  other  side  let  X  >  X^  (fig.  5).  The  normal  to  the 

line  X  =  X]^  is  drawn  in  the  direction  of  decreasing  velocity. 

By  considering  a  certain  point  of  this  line  and  by  taking  the 
direction  of  the  y  axis  to  agree  with  the  direction  of  the 
normal  to  this  point,  from  equation  (1.6) 

(1  -  Mj^2)  sin  e^^  1^  -  X^^  sin  ^^^  |^  +  \  cos  ^i  |^  =  0 

cos  01  1^  -  Xi  cos  01  ^  -  An  sin  9,  |^  =  0       (1.7) 
•^  dn    ^  ^  OS         ^  -'■on 

where  M]_  is  the  Mach  number  for  X  =  X]^,  and  9-^     is  the  angle 
between  the  velocity  vector  and  the  positive  direction  of  the 
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tangent  to  the  line  X  =  X^  at  a  fixed  point  of  this  line.  By- 
eliminating  from  equations  (l.7)  the  magnitude  bd/^n     and  by 
determining  from  the  obtained  equation  SQ/Ss 


M  _  1  -  V  s^n^  ^1  5X 
Ss         X,        Sn 


(1.8) 


2    2 
In  the  subsonic  flow,  1  -  M-^  sin  0^^  >  0  for  any  9     and 

because  ^X/hn  <  0,     then  S9/5s  <  0.  Thus  9,      in  this  case, 
changes  monotonically  along  the  line  X  =  constant. 

By  considering  the  transition  line  vhere  X,  =  1  and  by 

starting  from  the  assiunption  that  the  transition  line  is  not  a 
line  of  discontinuity  (reference  l),  that  is,  that  all  the  deriva- 
tives 90/Ss,  ^/Sn,  SX/Ss,  and  Sx/Sn  are  finite  on  the  tran- 
sition line,  the  following  equation  is  obtained  from  equation  (1.8) 

^  =  cos2  0  |Z^  (1.9) 

Ss        on 

and  therefore  the  condition  S©/Ss  <  0.  Thus 

THEOREM  2.  -  If  a  point  moves  along  the  transition  line  so 
that  the  region  of  subsonic  velocity  lies  to  the  left,  the  velocity 
vector  will  monotonically  turn  in  the  clockwise  direction. 

The  condition  B0/Ss  =  0  along  the  transition  line,  in  the 
case  where  a  transition  from  the  subsonic  to  the  supersonic  veloc- 
ity occurs,  was  previously  obtained  by  S.  A.  Christianovich. 

The  property  of  monotonicity  is  not,  however,  characteristic 
for  supersonic  flow.  In  this  case,  the  inequality  1  -  M^  sin^  0  >  0 
and  the  inequality  1  -  M^  sin^  0  <  0  could  hold  and  therefore  the 
value  of  S0/Ss  can  change  sign. 

From  further  discussion,  it  will  be  clear  that  the  fact 
expressed  by  theorem  2  determines,  to  a  considerable  extent,  the 
character  of  the  flow  in  the  local  supersonic  region  and  also  the 
possibility  or  impossibility  of  the  simultaneous  existence  of  the 
subsonic  and  supersonic  flows  without  change  in  the  potential  char- 
acter of  the  flow.  This  fact  will  hereinafter  be  termed  the  "law 
of  monotonicity  on  the  transition  line,"  or  simply  the  "monotonic 
law." 
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3.  Property  of  Monotonic  Change  of  Velocity  Magnitude 
and  Its  Inclination  along  Characteristics 

By  making  use  of  the  results  of  section  2,   it  will  be  proved 

that: 

THEOREM  3.  -  If  in  the  supersonic  region  there  is  given  a 
section  of  the  characteristic  of  one  family,  such  that  the  char- 
acteristics of  the  other  family  originating  from  the  points  of 
this  segment  end  on  the  transition  line,  the  angle  of  inclination 
of  the  velocity  vector  and  the  magnitude  of  the  velocity  are  mono- 
tonic  functions  along  the  given  segment  of  the  characteristic. 

By  considering  a  certain  point  C  on  the  segment  AB  of  the 
characteristic,  draw  from  the  points  A,  B,  and  C  to  the  line  of 
transition  the  characteristics  of  the  family  different  from  the 
family  to  which  the  cloaracteristic  AB  belongs  (fig.  6).  The  ends 
of  these  cliaracteristics  lying  on  the  transition  line  are  A^.,  B^, 
and  C^,  respectively. 

In  the  plane  of  the  hodograph,  points  A,  B,  and  C  are 
represented  by  the  points  A',  B',  and  C,  respectively,  lying 
on  a  certain  epicycloid  y;     and  the  points  A,^,  B^,  and  C^^ 

are  represented  by  the  points  A^ ' ,  B^ ' ,  and  C^ ' ,  respectively, 
lying  on  the  circle  X  =  1.  Each  pair  of  points.  A'  and  A,,.', 
B'  and  B^',  and  C  and  C^'   lie  on  one  of  the  epicycloids  of 

the  family  different  from  the  one  to  which  the  epicycloid  y 
belongs  (fig.  7). 

As  the  point  C  moves  along  the  epicycloid  y     and  its  cor- 
responding point  C^'  moves  along  the  circle  X  =  1,  the  polar 

angles  0  at  the  points  C '  and  C^ '   either  simultaneously 
decrease  or  simultaneously  increase  because  each  straight  line, 
9  =   constant,  intersects  the  epicycloid  only  at  one  point.  When 
point  C  moves  in  the  same  direction  from  point  A  to  point  B 
along  the  segment  of  the  characteristic  AB,  the  point  C^^  moves 
along  the  transition  line  from  point  A^  to  point  B^^  likewise 
in  one  direction  as  the  characteristics  of  one  family  do  not 
intersect.   In  accordance  with  theorem  2  of  the  preceding  section, 
it  follows  that  the  point  C#'   likewise  moves  along  the  circle 
A  =  1  in  the  same  direction.  Hence,  the  polar  angle  9     for  the 
point  C  varies  monotonically  and,  as  follows  from  the  properties 
of  the  epicycloid,  the  magnitude  of  the  velocity  X  likewise 
varies  monotonically  for  point  C. 
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In  the  supersonic  region  I  about  the  profile  determined  in 
section  1,   both  families  of  characteristics  satisfy  the  conditions 
of  theorem  3  so  that  in  moving  along  any  characteristic,  the  inclina- 
tion of  the  velocity  vector  0     and  the  magnitude  of  the  velocity  X 
vary  monotonically.  For  motion  directed  toward  the  transition  line 
along  the  characteristics  of  the  first  family,  both  magnitudes  9 
and  X  monotonically  decrease;  whereas  the  characteristics  of  the 
second  family  (the  left  one),  9     monotonically  increases  and  X 
monotonically  decreases.  From  the  obtained  property  of  the  mono- 
tonic  variation  of  the  magnitude  of  the  velocity  and  the  angle  of 
inclination  of  the  velocity  vector  along  each  of  the  character- 
istics of  region  I,  it  follows  that  this  region  is  represented  as 
a  single  sheet  on  the  corresponding  region  of  the  hodograph. 

In  the  supersonic  region  II,  only  the  characteristics  of  the 
first  family  satisfy  the  conditions  of  theorem  3.  On  moving  along 
the  characteristics  away  from  the  profile,  both  magnitudes  9     and 
X  monotonically  decrease.  On  moving  along  the  characteristics  of 
the  second  family  in  region  II,  the  property  of  monotonicity  of  the 
change  in  the  magnitudes  9     and  X  does  not,  in  general,  hold. 

By  making  use  of  the  known  properties  of  the  characteristics 
in  the  plane  of  the  flow  and  in  the  plane  of  the  velocity  hodo- 
graph, the  directions  of  concavity,  as  shown  in  figure  8,  are 
obtained  in  the  region  of  applicability  of  theorem  3  for  small 
Mach  niombers.  This  direction  of  concavity  of  the  characteristics, 
as  follows  from  the  results  of  S.  A.  Christianovich  (reference  1), 

takes  place  only  in  the  case  M  <  Mq  =  2/«/3  -  x  =  1.565.  By  mak- 
ing use  of  the  results  mentioned,  it  is  found  that  if  M  =  Mq  for 
a  certain  point  K  in  the  region  of  applicability  of  theorem  3, 
the  characteristics  have  the  appearance  represented  In  figure  9. 
Point  K  is  always  a  point  of  inflection  of  the  characteristics. 


4.  Variation  of  Velocity  at  Profile  in  Supersonic  Region 

It  is  assxjmed,  as  is  usually  the  case,  that  in  the  supersonic 
region  the  profile  is  everywhere  directed  convex  to  the  flow.   In 
the  supersonic  region  at  the  profile,  an  arbitrary  point  E  and 
the  characteristic  of  the  second  family  KKg*  are  assumed  to 
originate  at  this  point  (fig.  10).  For  point  K  from  equation  (l.l), 

is  the  function  introduced  according  to 


X     =  f(02» 

-^k) 

where 

f 

is  the  fu 

equation  (l 

.1). 

d( 

dX      _ 
-  ^k)   " 

f 

(^2*   -   \) 

^^2. 


'*d(T^ 


(1.10) 
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In  accordance  with  theorem  2,  6.62^/ i{-  \)  ^  0;     hence,  by- 
making  use  of  equation  (l.lO) 

^^^<f'(e2*  -0k)  (1-11) 

But  02*  "  ^k  ~  *P(^)^  where  cp  is  the  reciprocal  function 
of  f.  Hence 

f'(02*-ek)  =f'MA)3  =§ 

where  the  differentiation  is  effected  along  the  epicycloid  in  the 
plane  of  the  hodograph  (fig.  11).  Two  points  A  and  B  are 
assumed  to  he  on  the  epicycloid,  the  polar  angle  between  which  is 
equal  to  d0.  The  value  of  X  is  also  assumed  to  be  greater  at 
point  B  than  at  point  A. 

From  point  A,  an  arc  of  the  circle  X  =  constant  is  drawn  to 
the  intersection  with  the  straight  line  OB  at  point  C,  the 
length  of  the  arc  being  taken  equal  to  da.  The  angle  between  the 
tangents  to  the  arc  of  the  circle  AC  and  the  arc  of  the  epicy- 
cloid AB  at  point  A  is  equal  to  the  Mach  angle  a.  Hence, 
dX  =  tg  a  da  and  because  da  =  Xd9,  then  dV<i0  =  ^"tg  a. 

The  graph  of  the  function  f'[cp(X)]  =  Atg  a  is  given  in 
figure  5.  The  inequality  (l.ll)  therefore  assumes  the  form 

dFT7-  ^"^S  a        ■         (1.12) 

As  is  shown  by  the  preceding  inequality,  this  value  depends 
only  on  the  magnitude   X  at  the  given  point.  By  comparing  equa- 
tion (1.10)  emd  inequality  (l.ll),  it  is  found  that  inequality  (1.12) 
may  be  expressed  as  an  equation  if  and  only  if  at  a  given  point  of 
the  profile  dQg^/d^k  =  0*  I^  "this  relation  holds  true  on  a  cer- 
tain segment  AB  of  the  profile  (fig.  12),  the  condition  Qg*  = 

constant  is  obtained.  All  the  characteristics  of  the  second  fam- 
ily starting  from  AB  are  then  represented  in  the  plane  of  the 
hodograph  by  the  same  epicycloid  of  the  second  family  7   (fig.  13) 
because  at  their  ends,  lying  on  the  line  X  =  1,  the  velocity  is 
constant.  Hence,  the  entire  region  ABB„  A^    (fig.  12)   is  also 

represented  in  the  plane  of  the  hodograph  by  the  segment  of  the 
epicycloid  A'B',  upon  which  the  points  of  the  segment  AB  will 
also  lie. 
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The  characteristics  of  the  first  family  in  the  region  ABB„  A 

2*  o* 

must,  in  this  case,  be  straight  lines  with  constant  velocity  on 
each,  that  is,  in  this  region  a  certain  Meyer  flow  (rarefaction) 
must  occur.  However,  as  will  be  shown  in  section  6,  the  realiza- 
tion of  such  flow  under  the  conditions  of  the  problem  is  impossible. 
Hence,  in  the  relation  (1.12),  the  equality  can  hold  only  at  certain 
points  of  the  contour.  ~ 

In  the  supersonic  region  I,  the  characteristics  of  the  first 
family  likewise  end  on  the  line  X  =  1.  By  applying  considerations 
analogous  to  those  previously  mentioned,  the  second  inequality  along 
the  profile  in  region  I  is  obtained.- 

^^^>-Xtga  (1.13) 

In  the  case  where  the  equality  in  relation  (1.13)  is  attained 
on  a  certain  segment  of  the  profile  AB  in  the  region  bounded  by 
the  characteristics  of  the  first  family  originating  at  points  A 
and  B,   the  transition  line,  and  the  profile,  a  Meyer  flow  (compres- 
sion) teikes  place,  the  straight  lines  being  the  characteristics  of 
the  second  family  (fig.  14).  As  will  be  shown  in  section  6,  such 
flow  likewise  cannot  be  realized  under  the  conditions  of  the  problem 
and  therefore,  in  relation  (1.13),  the  equality  can  be  attained  only 
at  certain  points  of  the  contour.  By  combining  relations  (1.12) 
and  (1.13),  the  following  inequality  is  obtained  on  the  contour  of 
region  I:  • 


dX 


d(-  e.  ) 


<  Xtg  a  (1.14) 


Inequality  (1.13)  is  of  little  interest  in  the  case  of  a  shock 
wave  that  closes  the  supersonic  region  because  in  this  case  in 
region  I,  the  velocity  generally  does  not  decrease  so  that  the  rela- 
tion refers  mainly  to  the  case  of  the  flow  with  supersonic  region 
without  a  shock  wave.   In  this  case,  region  I  coincides  with  the 
entire  supersonic  region. 

The  equation  for  the  change  in  the  velocity  can  also  be 
obtained  for  the  case  where  a  certain  segment  of  the  profile  is 
concave  to  the  flow.  Thus,  for  any  point  C  on  the  segment  AB 
(fig.  15) 

^=f(^2*-^k)       i^  =  ^'(^2*-^k)(Sr-    ) 
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But  i"(^2-)(-  ~  ^k)  =  A  tg  a,  hence 

^<  -  Atga  (1.15) 

By  making  use  of  the  law  of  monotonia ity,  d©2^^/ci■^v  ^  0  and 
therefore  from  inequality  (1.15) 

^<  -  A  tga  •     (1.16) 

[NACA  Comment:  equation  (1.16)  is  apparently  in  error.] 

Thus,  on  the  segment  of  the  profile  having  a  concavity  facing 
the  flow,  the  velocity  drops  and,  as  follows  from  relation  (1.16), 
its  drop  cannot  be  too  slow.   Inequality  (1.16)  differs  essentially 
from  relations  (1.12)  and  (1.13),  which  show  that  on  the  convex 
segment  of  the  profile,  the  velocity  cannot  vary  too  rapidly. 

Strictly  speaking,  these  conditions  are  true  up  to  those 
changes  in  the  flow  that  are  brought  about  by  an  oblique  shock 
wave  (fig.  15).  By  considering  "the  characteristic  of  the  second 
family  ^2*^1      ^^^   ^^^  prolongation   C]_C,   it  is  found  that  the 
transition  through  the  shock  wave  in  the  flow  plane  corresponds 
to  the  displacement  along  the  segment  of  the  strophoid  C'-,C., 

in  the  plane  of  the  velocity  hodograph  (fig.  16).  Hence,  the 
transforms  of  the  segments   C^Cp    and   CnC   of  the  character- 
istic lie  on  two  different  epicycloids  of  the  second  family.  How- 
ever, by  making  use  of  the  fact  that  at  the  point  C'l   the  stro- 
phoid and  epicycloid  have  a  common  tangent  and  the  same  radius  of 
curvature,  it  is  found  that  with  an  accuracy  up  to  small  magni- 
tudes of  the  third  order,  relative  to  the  changes  in  velocity  in 
the  shock  wave,  these  two  epicycloids  may  be  considered  as  coin- 
ciding. Inequality  (1.16)  holds  with  the  same  degree  of  accuracy. 


5.  Flow  in  Supersonic  Eegion  Arising  from  Presence 

of  Straight  Segment  on  Profile 

The  straight  segment  AB  is  assiuned  to  be  on  the  profile  in 
the  supersonic  region.  The  characteristic  of  the  second  family 
CCg*  originates  from  the  arbitrary  point  C,  which  is  assumed  to 

be  on  the  straight  segment  AB  (fig.  17).  At  point  C,  the 
equation  is  X=  f(02«  "  ®k)* 
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As  point  C  moves  along  segment  AB  in  the  flow  direction, 
the  magnitude  6,      remains  constant  because  of  the  rectilinearity 

of  the  segment  and  the  magnitude  ^g*'   ^^  accordance  with  theo- 
rem 2,  does  not  increase.  Hence,  using  the  fact  that  f  is  a 
monotonically  increasing  function  of  its  argument: 

THEOREM  4 .  -  On  a  straight  segment  of  a  profile  in  a  super- 
sonic region,  the  velocity  in  the  flow  direction  does  not  increase. 

It  is  assumed  that  the  straight  segment  is  located  in  the 
supersonic  region  I  (fig.  18). 

From  the  ends  of  the  straight  segment  to  the  transition  line, 
the  characteristics  of  the  first  family  fJ^i*     and  BBj^^  and  of 

the  second  family  AA2^(.  and  BB2^  are  drawn.  Let  9-^  =  9q     on 

AB.  From  point  C  near  point  A,  the  characteristic  of  the 
second  family  is  drawn  to  the  intersection  with  the  character- 
1  stic  AA]_^  at  point  C.  From  theorem  3,  it  is  found  that 

e(A)  >0(C')  >0(C),  but  0(A)  =  0(C)  =  00,  hence  0(C')  =  0o 
and,  as  again  follows  from  theorem  3,  0  =  0o  over  the  entire 
segment  of  the  characteristic  CC. 

By  moving  point  C  from  point  A  to  point  B,   it  is  found 
that  the  characteristics  AAy^.     and  BB2#  necessarily  intersect 
and  in  the  triangle  ABD,  formed  by  the  segment  AB  and  these 
characteristics,  0  =  0o  =  constant.  Hence,  it  follows  that 
X  =  constant  and  the  characteristics  of  both  families  are 
rectilinear. 

The  entire  region  ABBp  A_    is  represented  in  the  plane  of 
the  velocity  hodograph  by  a  single  epicycloid  of  the  second  family, 
as  all  the  characteristics  of  the  second  family  are  represented  by 
a  single  epicycloid  passing  through  the  point  0  =  6q,      X  =  A(A). 

Hence,  over  the  entire  region  considered,  a  certain  Meyer  flow 
(rarefaction)  will  take  place  with  straight  characteristics  of  the 
first  family.  Similarly,  the  entire  region  ABB,  A^    is  repre- 
sented in  the  plane  of  the  hodograph  by  a  single  epicycloid  of  the 
first  family,  over  which  a  Meyer  flow  (compression)  will  take 
place  with  straight  characteristics  of  the  second  family. 

As  was  previously  stated  and  as  will  be  shown  in  section  6, 
the  realization  of  Meyer  flows  under  the  conditions  of  the  problem 
is  lii5)ossible. 
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6.  In^jossibility  of  Meyer  Flow  in  Local  Supersonic  Region 

In  the  plane  x,y,  a  Meyer  flow  is  assxjmed  to  "be  between  the 
characteristics  C-|_  and  C2  of  the  second  family  with  straight 

characteristics  of  the  first  family  (fig.  19). 

From  a  certain  point  A  of  characteristic  Ci,     the  straight 
characteristic  of  the  second  family  is  drawn  to  its  intersection 
at  point  B  with  the  characteristic  Cg.  On  the  straight  line  AB, 

the  Mach  angle  a  =  constant,  hence,  it  intersects  the  character- 
istics Ci  and  Cg  at  the  same  angle  7  =  «  -  2a, 

Set  AB  =  Z  and  denote  by  3  the  angle  formed  by  the  straight 
line  AB  with  the  ax3.s  of  abscissas.  The  straight  characteris- 
tic A]^B]_,   infinitely  near  AB,   is  considered.  The  segments  of 

the  straight  lines  AB  and  AiB]_  are  denoted,  respectively,  by 

dS]_   and  dsg   on  the  characteristics  C^^  and  Cg  and  the  angle 

between  AB  and  A^Bi  is  denoted  by  d3.  Thus,  the  following 
equations  are  obtained: 

ds,  =  ^1^     ds2  =  Ll-LllM  (1.17) 

1   sin  7       ^  sin  7  ^    ' 

where  r  is  the  distance  between  point  A  and  the  point  of  inter- 
section of  the  straight  lines  AB  and  AiBi. 

Set  the  length  of  the  segment  AiBi  equal  to-  I  +  dZ .  From 
figure  18,  the  equation  dZ  =  (dsg  -  dsi)  cos  7  is  obtained,  or 
by  making  use  of  equations  (1.17) 

dZ  =  Z  ctg  7  dp  (1.18) 

In  the  plane  of  the  hodograph,  the  characteristic  Ci  is 
represented  by  an  epicycloid. 

The  directions  AB  and  AiBi  coincide  with  the  direction  of 
the  normals  to  this  epicycloid  drawn  at  the  points  A'  and  A'l, 
corresponding  to  the  points  A  and  Ai  (fig.  20).  The  angle 
between  the  radius  vector  and  the  direction  of  this  normal  is  equal 
to  the  Mach  angle  a.  Hence,  p  =  a  +  9  and  dp  =  da  +  d0. 

Thus,  by  making  use  of  equation  (1.18) 

^  =  ctg  7 (da  +  de)  =  -  ctg  2a  da  +  ctg  7  d0      (1.19) 


14  NACA  TM  No.  1213 

or  "by  integrating 

vhere  Iq,     clq,     and  9q  characterize  a  certain  initial  point. 

For  a  >  «/4  and  9  >  9q,     both  terms  on  the  right  of  equa- 
tion (1.20)  are  positive  so  that  the  following  inequalities  hold: 


m  f  >  in  /i^^  I  >  2o  /^^^  (1.21) 

Iq  V  am  2a  ^  ^  sin  2a 

If  X  ^  1,  then  a  -^  rt/2  and  Z  ->». 

An  entirely  analogous  result  is  obtained  for  the  Meyer  flov 
betveen  the  characteristics  of  the  first  family  with  straight 
characteristics  of  the  second  family.  Thus 

THEOEEM  5 .  -  A  Meyer  flow,  between  two  characteristics  of  one 
family  with  straight  characteristics  of  the  second  family,  cannot 
be  entirely  realized  up  to  the  line  X  =  1  in  a  finite  region. 

If  on  the  characteristic  C-^,     represented  in  figure  19,  the 

equality  A  =  1  is  attained  at  a  certain  point  D,  then  by  making 

use  of  the  preceding  result,  it  is  found  that  the  tangent  at 

point  D  to  the  characteristic  C-,      is  asymptotic  for  any  char- 

■^  1 

acteristic  C2  placed  relatively  to  C]_,  as  shown  in  figure  19  . 

From  theorem  5,  it  follows  that  along  the  straight  part  of  the 
profile  in  the  local  supersonic  region,  the  magnitude  of  the  veloc- 
ity cannot  be  constant.  If  such  were  the  case,  then  between  the 
characteristics  of  the  second  family  originating  at  the  ends  of 
this  segment  a  Meyer  flow  would  take  place.  This  realization, 
because  of  the  finiteness  of  the  local  supersonic  region,  is 
impossible  as  follows  from  theorem  5.  Thus  theorem  4  may  be  more 
definitely  stated  as  follows: 

^According  to  S.  A.  Christianovich,  if  a  Meyer  flow  originates 
from  a  straight  characteristic  on  which   A  =  1,  the  cliaracter- 
istics  of  the  other  family  do  not  originate  from  this  straight 
line. 
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THEOEEM  6.  -  On  the  straight  segment  of  a  profile  in  the  super- 
sonic region^  the  velocity  decreases. 

II  -  CONDITIONS  FOE  BREAKDOWN  OF  POTENTIAL  FLOW 

WITH  LOCAL  SUPERSONIC  REGION 

The  occiirrence  of  a  wave  resistance  for  a  flov  with  large  sub- 
sonic velocity  about  a  body  is  assumed  to  be  connected  with  the 
instant  when  at  any  point  on  the  surface  of  the  body  a  velocity  is 
reached  equal  to  the  local  velocity  of  sound. 

The  study  of  the  results  of  a  number  of  tests  has  shown,  how- 
ever, that  the  occurrence  of  a  shock  wave  and  therefore  the  arising 
of  a  wave  resistance  sometimes  takes  place  beyond  Mcr*  Photographs 
obtained  by  the  Toepler  method  show  In  certain  cases  a  local  super- 
sonic region  extended  deep  into  the  flow  without  any  shock  wave. 

The  necessity  for  the  breakdown  of  the  potential  character  of 
the  flow  with  local  supersonic  region,  that  is,  the  necessity  for 
the  appearance  of  a  shock  wave,  has  not  yet  been  theoretically 
demonstrated  for  a  single  case.  It  will  therefore  be  of  interest 
to  describe  certain  conditions  that  are  known  to  be  accompanied  by 
a  breakdown  of  the  potentiail  flow, 

1.  Deformation  of  Contour 

-  How  the  deformation  of  the  contour  can  lead  to  the  breakdown 
of  a  previously  existing  potential  flow  is  considered. 

A  flow  with  local  supersonic  region  about  a  curvilinear  con- 
tour for  a  Mach  number  of  the  oncoming  flow  such  that  the  flow  is 
potential  (fig.  21)  is  considered. 

From  points  A  and  B  of  the  contour,  the  characteristics 
of  both  families  are  drawn  both  in  the  flow  plane  and  in  the  plane 
of  the  hodograph.  Because  of  satisfying  the  condition  of  monotonic- 
ity  of  change  of  the  angle  on  the  transition  line,  the  transforma- 
tions of  the  points  Ag^,  Bg*,  f^i^f     and  B^^^.  in  the  hodograph 

plane  axe  located  on  the  circle  A  =  1  in  the  same  order  as  on  the 
transition  line  in  the  flow  plane.  The  contour  is  deformed  in  such 
a  way  that  the  segment  of  the  arc  AB  transforms  into  a  straight 
segment  (fig.  22). 
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By  virtue  of  theorem  6,  section  1,   the  velocity  along  the 
straight  segment  AB  can  only  decrease.  By  drawing  the  character- 
istics in  the  hodograph  plane  and  by  making  use  of  theorem  1,  it 
is  found  that  the  transforms  of  the  points  Ag^^,  ^2*>     ^1*'  ^^^ 
Bi»  are  not  located  on  the  circle  X  =  1  in  the  order  in  which 
these  points  are  arranged  on  the  transition  line  in  the  flow  plane. 
Consequently,  the  deformation  of  the  contour  breaks  up  the  mono- 
tonicity  on  the  transition  line.  Thus 

THEOEEM  7.  -  The  characteristics  of  the  first  family  that 
originate  from  the  points  of  a  straight  segment  of  the  contour  can 
never  end  on  the  transition  line  but  must  fall  on  a  shock  wave. 

If  the  contour  deformation  considered  is  effected  continuously, 
it  is  evident  that  the  breakdown  of  the  potential  flow  occurs  before 
the  deformed  segment  becomes  straight.  This  fact  permits  the  con- 
clusion that  at  a  given  Mach  niimber  of  the  oncoming  flow,  any  pro- 
file convex  to  the  flow  can  be  so  deformed  that  a  new  profile  is 
again  obtained  convex  to  the  flow,  the  flow  about  which  for  the 
same  Mach  number  is  no  longer  possible  without  a  shock  wave. 


2.  Impossibility  of  Formation  of  Line  of  Discontinuity 

Within  Local  Supersonic  Region 

Investigation  will  reveal  whether  the  impossibility  of  realiza- 
tion of  a  continuous  flow  within  the  supersonic  region  for  given 
boundary  conditions,  as  expressed  in  the  intersection  of  the  char- 
acteristics of  one  family  and  the  formation  of  a  supersonic  shock 
wave  within  the  region,  can  be  a  cause  for  the  breakdown  of  the 
potential  character  of  the  gas  motion. 

Karman  (reference  3)  makes  the  assumption  that  the  probable 
cause  of  the  appearance  of  shock  waves  is  the  formation  of  an 
envelope  of  the  Mach  lines  in  the  supersonic  zone  for  a  certain 
Mach  mrniber  of  the  oncoming  flow. 

The  envelope  of  the  characteristics  of  one  family,  as  shown 
by  S.  A.  Christieuaovich  (reference  l),  coincides  with  the  line 
of  discontinuity  that  is  determined  as  the  geometric  locus  of  the 
points  at  which  at  least  one  of  the  derivatives  S9/Sx,  Ss/Sy, 
h}\/^x,     or  ^V^y  becomes  infinite.  Hence,  to  investigate  the 
problem  of  the  possibility  of  the  formation  of  supersonic  shock 
waves  within  the  supersonic  region  on  increasing  the  Mach  number 
of  the  oncoming  flow,  it  is  necessary  to  obtain  the  value  of  the 
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acceleration  at  the  points  of  the  supersonic  region.  It  will  be 
shown  that  on  the  contour  of  the  profile  infinite  accelerations 
ceuinot  occur. 

Given  a  certain  flow  in  the  supersonic  region  with  continuous 
transition  from  supersonic  to  subsonic  velocities  (sec.  I,  pt.  1, 
case  (a));  assume  a  certain  point  K  on  the  contour  in  the  super- 
sonic region  and  take  the  x  axis  in  the  flow  direction  (fig.  23). 
Eqviatlons  (1.6)  then  become 


5X 


(m2-i)^-x|£  =  o 


(2.1) 


Thus 


BX   dX  Se 


^ 


^a^  Sx 


dX 


k  =  -  i 


1  dX_ 
E  6.9, 


(2.2) 
(2.3) 


where  k  is  the  curvature  of  the  contour  at  point  K  and  E 
the  radius  of  the  curvature  at  this  point.  By  making  use  of 
Inequality  (1.14)  and  equations  (2.2)  and  (2.3) 


is 


ax 


^  I  >v  tg  a 


he 


1 

E 


(2.4) 


From  equations   (l.20) 


Be 


<  ^  ctg  a 


ax 


E 


(2.5) 


Eelations  (2.4)  and  (2.5)  show  that  the  magnitudes  hd/hx, 
hO/bj,     h'h/bx,     and  BV^7  ^^e  always  finite  on  the  contour  within 
the  supersonic  region  provided  that  everywhere  E  ^  0.  The  case 
E  =  0  (infinite  curvature)  at  some  point  is  a  special  case  and  never 
occurs  on  real  profiles.  Therefore,  assxime  that  on  the  contovirs 
considered,  at  all  points  E  >  Eq  where  Eq  is  a  certain  constant. 

The  maximum  possible  value  of  the  acceleration  at  the  interior 
points  of  a  local  supersonic  region  is  computed.  The  differential 
equations  of  the  characteristics  of  the  first  and  second  families 
in  the  plane   x,y   are,  respectively. 
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dy  =  tg(e  +  a)  dx       dy  =  tg(0  -  a)  dx       (2.6) 


Moreover,  along  the  characteristics  of  the  first  and  second 
families,  respectively,  the  relations  hold  (see  for  example 
reference  1) 


a  -  0  =  2ti  =  constant    (2.7) 


^^  -  ^         4^      (2.8) 


C(x-  i)/(x+  1)3x2  X 


From  relations  (2.6)  and  (2.7) 

By  differentiating  the  first  relation  with  respect  to  |,  the 
second  with  respect  to  t),  and  by  subtracting  one  from  the  other, 
the  following  equation  is  obtained  for  the  value  of  x: 

JL ii|^  |x  ^g_^j  1^  ^  _i ^ie^  |x  ^  ^g(0^)  ^ 

co82(0-a)       ^Ti       5|         ^'        '  SlSn       co82(0m)       a  I       Sti  SISti 

(2.10) 

From  relations   (2.7) 

a    =    |+Ti  e=|-Ti  (2.11) 

By  making  use  of  these  equations  and  choosing  the  coordinate 
system  so  that  the  direction  of  the  x  eixis  coincides  with  the 
direction  of  the  velocity  at  some  point  A  (fig.  24),  the  relation 
between  the  partial  derivatives  of  x  at  the  point  considered  is 
obtained  from  equation  (2.10) 

S^x   1  +  a'  (C)  hx      1  +   a'(a)  5x  _        /  t(^\.^ 
^I^T    sin  2a  hi         sin  2a   ^  "        ^'^  ^^'   ~  do 

(2.12) 
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Let  dZ  and  ds  "be  the  lengths  of  the  elementary  segments  of 
the  characteristics  of  the  first  and  second  families,  respectively. 
In  passing  from  point  A  to  neighboring  points  in  the  directions 
indicated  in  figure  24,  dl  and  ds  will  be  taken  as  positive 


Sx  dl 
hi   d| 


cos  a 


dZ 
dl 


Bx  Sx  ds  ds  /„  ,_v 
•^-  =  •^-  -T—  =  cos  a  -T—  (2.13) 
dTJ   OS  dii        dT]  ^    ' 


vhere  the  derivative  dl/dt,    is  taken  along  a  characteristic  of  the 
second  family  and  the  derivative  ds/dri  is  taken  along  a  char- 
acteristic of  the  first  family.  By  making  use  of  equations  (2.13), 
relation  (2.12)  can  be  reduced  to  one  of  the  two  following  forms: 


^  r "  °'  df ) 


^1 


cos  a 


1  +  g' (g) 
sin  2a 


1  +  a'(a) 
sin  2a 


cos  a 


dl^ 


cos  a 


dTl; 


1  -H  a'(a)  ds 
2  sin  a  dT) 


1  + 


a'(a)  dl 
2  sin  a  d| 


=  0 


=  0 


(2.14) 


(2.15) 


Along  the  characteristic  I   =  constant,  relation  (2.14)  may 
be  considered  as  an  ordinary  differential  eqviation  relative  to  the 
magnitude  cos  a  dZ/d|  and  similarly  along  the  characteristic 
Tj  =  constant,  the  relation  (2.15)  may  be  considered  as  an  ordi- 
nary differential  equation  relative  to  cos  a  ds/dri.  By  setting 
z  =  cos  a  dZ/d^,  the  following  linear  equation  of  the  first  order, 
based  on  equation  (2.14),  is  obtained  for  determining  z  along 
the  characteristic  i,   =  constant: 

dz  ^  l-Hx'(a)  ^  .  l^a'(a)  dz  ,    . 

dTi   sin  2a     2  sin  a  dii  K^.i-b) 

The  solution  of  this  differential  equation  has  the  form 


z  =  exp 


/ 


L  ''0 


^4^dri 
sin  2a   ' 


/ 


exp 


/ 


1+a'  (g) 
sin  2a 


dTl 


l+g' (g)  ds  ,  , 
2  sin  a  d^  ^ 

(2.17) 


where  t^q  is  the  value  of  t\     at  a  certain  initial  point  on  the 

characteristic  considered,   |  =  constant,  and  Zq  =  cos  aQ(dZ/d|)o 

is  the  value  of  z  at  this  point.  The  integral  under  the  sign  of 
the  exponential  function  in  equation  (2.17)  can  be  readily  computed. 
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if  the  magnitudes  entering  the  expression  under  the  integral  are 
expressed  in  terms  of  the  Mach  angle  a.,     hy  making  use  of  the 
relation 


X  = 


..E 


ysin^  a  +  (x-l)/2 


and  the  relations 


dT]  =  da  =  -  de 


^  =  X  tg  a 


which  hold  along  the  characteristic  of  the  first  family  I,  =  constant, 
the  first  of  which  follows  from  relations  (2.7)  and  (2.11)  and  the 
second  of  which  was  already  obtained  in  section  I,  part  4.  The 
following  expression  is  thus  obtained: 


/ 

^0 


sm  2a    '      I  V'tg  CO, 


'0 


6in2^  ^Ir?   Oq  +  1/2  (x-1)] 


sin^ttQ  [sin^  a  +  1/2(^-1)] 


(2.18) 


Hence,  by  replacing  in  eqviation  (2.16)  the  magnitude  z  by  its 
expression,  and  denoting  by  G(a,aQ)  the  expression  under  the 

logEorlthm  sign  in  (2.18) 


cos  a 


dl 


df   G(a,aQ) 


cos 


°X)  (H)  ^   J   G(a,ao) 


l-Hx'(a)  ds 
2  sin  a  dT) 


dTl 


(2.19) 


An  entirely  similar  solution  of  equation  (2.15)  gives 


cos  a 


ds 


dT)   G(a,aQ) 


/ds\     r   _/    V  l+a'  (o)   dl  ^  . 


'0   To 


(2.20) 


By  making  use  of  the  property  of  mono tonicity  of  change  of  the 
magnitude  of  the  velocity  and  the  angle  of  inclination  of  the 
velocity  vector  along  the  characteristics  in  the  local  supersonic 
region  (sec.  I,  pt.  3),  it  is  noted  that  dl/d^  >  0  and  ds/dT]  >  0. 
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Moreover,  a' (a)  >  0.  Hence,  from  equations  (2.19)  and  (2.20), 
the  following  Ineqxialities  are  obtained: 


dl 


cos  a  ^  >  cos  a^  (j^J^  -^^^-^ 


dl\ 


(2.21) 


ds  ^ 
cos  a  -r-  >  cos  a-, 

dT)       0 


ds 
dii 


1 


(2.22) 


In  order  to  obtain  the  accelerations  at  a  certain  point  K  of 
the  local  supersonic  region,  the  characteristics  of  the  first  and 
second  families  are  dravn  from  this  point  up  to  the  intersection 
with  the  contour  at  points  A  and  B,  respectively  (fig.  25). 
Let  a  =  a_^^   and  a  =  a„p  at  points  A  and  B,  respectively. 


01 


"02 


Also,  let  the  radii  of  curvature  at 
respectively. 


A  and  B  be  E,   and  E 


2' 


The  ineqiialities  (2.21)  and  (2.22)  at  the  point  K  are  con- 
sidered where  point  B  is  taJcen  for  the  initial  point  in  inequal- 
ity (2.21)  and  in  inequality  (2.22)  point  A  is  used.  By  making 
use   of  relations  (1.14),  (2.1),  and  (2.7),  approximations  are 
easily  obtained  for  the  magnitudes  d  i/dl     at  point  B  and  dr\/ds 
at  point  A 


§1 
dl 


B 


Ec 


cos  a. 


•02 


dn 


ds 


cos 


"Ol 


(2.23) 


By  taking  the  direction  of  the  x  axis  at  point 
streamline,  the  following  equations  axe  obtained: 


K  along  a 


he 

^=  - 


2  cos  a 


dZ   ds 


^  _  _  Xtg  g 
o^     2  cos  a 


'da 

ds 


61 


he 

2 

1 

sin  a 

\dl 

dsj 

SA 

- 

A  tg  a 
2  sin 

a  V <i^ 

ds 

(2.24) 


where  the  relation  a'(X)  =  -  ctg  a/ A  is  used.  By  malcing  use  of 
inequalities  (2.21),  (2.22),  and  (2.23)  and  equations  (2.24)  at 
point  K,  the  following  inequalities  holding  at  this  point  are 
finally  obtained: 
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i|  ^  ^^%1'\^   ^  ^^^02'^2[]  ^(^^  -^^^ 
1^1  <  [L(aQi,E3^)  +  1(002,^2^  N(a)7ctg  ( 


<  IlCoq^^R-]^)  +  L(aQ2*^2)l  ^^°-^   ^  *6  a  -/tg  a 


1^  <  fL(aQi,E-L)  +  L(aQ2,E2)l  N(a)X  TtilT      (2.25) 


where 


N(a)  = 


.  2 

sm  a 


jin'^a  +  1/2  (x-1). 


1 
2(X-1) 


L(aQ,E)  = 


NCoq) 


E 


•/ctg 


°0 


The  ineqxialities  (2.25)  shov  that  the  absolute  values' of  the 
derivatives  cannot  exceed  certain  finite  values,  depending  on  the 
Mach  number  at  the  given  point  of  the  supersonic  region  eind  the 
minimum  fixed  radius  of  curvature  on  the  part  of  the  contour  that 
is  a  boundary  of  the  supersonic  region. 

The  assumption  that  the  cause  of  the  formation  of  a  shock  wave 
is  the  envelope  of  the  characteristics  within  the  local  supersonic 
region  is  therefore  invalid.  Thus,  it  is  possible  to  fonaulate 

THEOEEM  8.  -  If  for  certain  conditions  determining  the  motion 
of  the  gas^  a  potential  flow  existed  with  local  supersonic  region, 
the  formation  of  a  shock  wave  within  or  on  the  boundary  of  the 
supersonic  region,  arising  from  any  change  in  these  conditions, 
cannot  be  preceded  by  the  occurrence  of  infinite  accelerations  at 
the  interior  points  of  the  supersonic  region.  That  is,  if  the 
curvatiire  of  the  part  of  the  contour  lying  in  the  supersonic  region 
does  not  become  infinite. 


'^When  the  flow  about  the  body  is  infinite,  the  conditions 
determining  the  motion  of  the  gas  are  the  Mach  number  of  the  oncom- 
ing flow  and  the  shape  of  contour  of  the  body. 
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The  theorem  obtained  refers  only  to  the  internal  points  of  the 
supersonic  region.  From  ineqiialities  (2.25),  no  conclusion  can  be 
drawn  as  to  the  finiteness  of  the  accelerations  on  the  transition 
line  because  the  right-hand  sides  of  three  of  them  approach  infin- 
ity as  the  transition  line  is  approached.  Hence,  the  question 
whether  the  shock  waves  are  preceded  by  the  occurrence  of  infinite 
accelerations  at  the  points  of  the  transition  line  remains  open. 


3.  Criterion  of  Collapse  of  Potential  Flow 

In  section  II,  part  1,  it  was  shown  that  a  deformation  of  the 
contour  can  lead  to  the  collapse  of  the  potential  flow.  Before  the 
deformation  of  the  curvilinear  segment  of  the  contour  at  its  points, 
the  inequality  d.}\/d.{-6-^)  >   -  Xtg  a  holds.  At  the  end  of  the 

deformation,  however,  at  the  points  of  the  obtained  straight  seg- 
ment, d.}\/d{-6-^)   =  -  00,  as  follows  from  theorem  6. 

Hence,  for  any  intermediate  state  of  the  deformed  contour 
there  is  first  attained,  at  some  point  of  the  segment,  the  equation 

^(^g  )  =  -  ^tg  a  (2.26) 

This  state  of  the  deformed  contour  is,  in  a  certain  sense, 
critical  because  for  further  defonaation  of  the  contour,  inequal- 
ity (1.13)  breaks  down  and  thus  the  flow  becomes  impossible  with- 
out a  shock  wave. 

In  the  given  case  condition,  equation  (2.26)  is  the  criterion 
of  breakdown  of  the  potential  flow.  Equation  (2.26)  may  also  be  a 
criterion  for  the  breakdown  of  the  potential  flow  in  the  case  where 
the  flow  is  about  a  fixed  contour  but  where  the  Mach  number  of  the 
oncoming  flow  increases. 

If  for  a  certain  Mach  number,  a  potential  flow  existed  with 
local  supersonic  region  and  at  the  points  of  the  contour,  the 
inequality 

<^^   >  -  X  tg  a 


d(-ej,) 


were  satisfied,  then  if  with  increase  in  the  Mach  number  of  the 
oncoming  flow  inequality  (1.13)  starting  from  a  certain  Mach  number 
breaks  down,  this  limiting  Mach  number  would  correspond  to  the 
attaining  of  equation  (2.26)  at  some  point  of  the  conto\ir. 
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The  Mach  number  of  the  oncoming  flow,  for  which  at  some  point 
of  the  contour  condition  (2.26)  is  satisfied,  will  be  denoted  as 
the  breakdown  Mach  number  M*-. 

The  nvunber  M*  is  thus  the  limiting  Mach  number  beyond  which 
the  body  is  subject  to  wave  resistance. 

The  breakdown  criterion  (2.26)  can  be  conveniently  repre- 
sented in  another  form,  by  considering  the  magnitude 

where  cp(X)  is  a  known  function  shown  in  figure  4.  By  differen- 
tiating this  equation  along  the  contour 


ds 


d0k 
ds 


But  cp'(X,  )  =  l/(X,  tg  cu  )  so  that  the  following  equation  for 
the  breakdown  criterion  is  obtained  equivalent  to  (2.26): 

de 

-^   =  0  (2.28) 

It  should  be  noted  that  the  breakdown  of  inequality  (1.13)  at 
any  point  of  the  contoiir  does  not  lead  to  the  breakdown  of  the 
potential  cheiracter  of  the  flow  near  this  point  but  makes  impossible 
the  ending  of  the  characteristics  of  the  first  family  from  the 
neighborhood  of  this  point  on  the  transition  line  as  the  law  of 
monotonicity  would  otherwise  fail  to  hold  at  the  transition  line. 
The  ch^aracteristics  of  the  first  family  must,  in  this  case,  merge 
in  the  shock  wave. 

In  order  to  determine  the  magnitude  M*  of  a  given  contour, 
it  is  necessary  to  know  how  for  an  increase  in  the  Mach  number  of 
the  oncoming  flow,  the  value  of  the  velocity  at  the  points  of  the 
contour  changes  for  a  potential  flow  with  local  supersonic  region. 

At  the  present  time,  certain  methods  are  known  for  the  approxi- 
mate solution  of  the  problem  of  the  flow  about  a  fixed  contour  for 
Mach  numbers  exceeding  M^j..   (See,  for  example,  references  4  and  5.) 
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In  all  these  methods,  the  convergence  of  the  process  of  the 
successive  approximations  has  not  been  proven;  hence,  the  results 
of  the  computations  made  with  the  aid  of  these  methods  must  be 
regarded  with  reserve.  It  seems  probable,  however,  that  up  to 
the  time  when  the  potential  flow  is  actually  possible,  these  methods 
at  least  qualitatively  represent  the  true  character  of  the  change 
in  the  distribution  of  the  velocity  on  a  fixed  contour  with  increase 
in  the  Mach  number  of  the  oncoming  flow. 

It  is  thus  of  interest  to  see  whether  from  these  solutions 
there  is  a  tendency  toward  the  attainment  on  the  contour  of  the 
condition  (2.28)  with  increase  in  Mach  number. 

From  the  geometric  data  of  a  wall  and  distribution  of  the 
velocity  along  the  wall,  obtained  in  the  work  of  Gortler  (refer- 
ence 5)  for  a  Mach  number  of  the  oncoming  flow  M  =  0.9,  the  angles 
of  inclination  0,^  were  computed  using  equation  (2.27).  The 

pattern  of  the  lines  of  flow  and  the  lines  of  equal  velocities  for 
this  case  is  shown  in  figure  26.  The  dependence  of  O^^     on  the 

coordinate  along  the  wall   |  is  represented  in  figure  27. 

As  may  be  seen,  the  flow  is  near  the  breaicdown.  The  monotonic 
decrease  of  the  angle  proceeds  up  to  i  =   5.5,  after  which  the 
angle  remains  almost  constant  or  slightly  decreasing. 

In  figure  28  is  shown  the  dependence  of  9^^     on  x  for  the 

flow  about  the  profile,  which  is  considered  in  an  American  inves- 
tigation as  M  =  0.75  and  M  =  0.83.  The  velocity  distribution 
over  the  contour  and  the  boundaries  of  the  local  supersonic  region 
are  shown  in  figure  29.   In  figure  28,  it  is  very  clearly  seen 
that  the  tendency''  toward  a  breakdown  of  the  monotonic ity  on  the 
transition  line  increases  with  increase  in  the  velocity.  At 
M  =  0.75,   somewhat  exceeding  M^j.,  the  angles  decrease  almost 
as  fast  as  on  the  contour.  For  M  =  0.83,  the  monotonicity  of 
change  of  the  angle  6-,        has  already  broken  down.  Hence,  at 

M  <  0.83  there  is  already  a  breakdown  of  the  potential  flow  near 
this  profile.  Strictly  speaking,  the  breakdown  for  M  =  0.83  no 
longer  has  significance.  Thus,  in  the  given  solutions,  a  tendency 
is  revealed  to  attain  on  the  contour  condition  (2.28)  on  increas- 
ing the  Mach  nvunber  of  the  oncoming  flow. 
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4.  Infinite  Accelerations  in  Accurate  Solutions 

In  section  II,  part  2,  it  is  shovn  that  an  infinite  accelera- 
tion on  a  contour  can  occur  only  if  the  curvature  of  the  contour 
at  any  point  is  infinite. 

An  exact  solution  of  a  certain  flow  with  local  supersonic  region 
with  the  occurrence  of  infinite  acceleration  at  two  symmetric  points 
of  the  contour  was  given  by  Eingleb  (reference  6).  According  to  the 
results  of  this  section,  this  flow  indicates  an  infinite  curvature 
of  the  contoiir  at  these  points.  Figure  30  shows  the  flow  pattern 
and  the  velocity  distribution.  At  points  P  and  Q  of  the  contour, 
Ringleb  obtained  infinite  values  of  the  derivative  of  the  velocity, 
that  is,  infinite  accelerations.  It  is  found  that  in  the  plane  of 
the  hodograph,  at  the  points  corresponding  to  P  and  Q  in  the 
flow  plane,  the  transform  of  the  streamline  and  the  epicycloid  have 
a  common  tangent. 

In  the  work  of  Karman  during  1941  (reference  3),  the  Eingleb 
solution  is  analyzed  in  detail  with  the  object  of  showing  for  this 
type  of  flow  the  possible  reasons  for  the  formation  of  -shock  waves 
in  the  local  supersonic  regions.  Karman  showed  that  if  at  the  con- 
tour in  the  supersonic  region  there  is  a  point  with  infinite  accel- 
eration and  a  value  of  S\i'/B9  f   w  (where  \j/  is  the  stream  func- 
tion and  0  is  the  angle  of  inclination  of  the  velocity  vector) 
then,  as  in  the  case  of  Eingleb,  the  transform  of  the  streamline 
in  the  plane  of  the  hodograph  touches  the  epicycloid  at  a  point 
corresponding  to  the  point  with  infinite  acceleration  in  the  flow 
plane.  Karman  identifies  this  contact  in  the  hodograph  plane  with 
the  presence  of  an  infinite  acceleration  at  the  corresponding  point 
of  the  flow  plane.  He  does  not,  however,  observe  the  fact  previously 
pointed  out  that  the  presence  of  an  infinite  acceleration  at  any 
point  of  the  supersonic  flow  indicates  the  presence  in  it  of  an 
infinite  curvature  of  the  streamline,  and  associates  the  infinite 
acceleration,  in  this  particular  case,  with  the  impossibility  of  a 
continuous  flow  in  the  general  case.  This  lack  of  observation  is 
indicated  by  the  fact  that  with  the  aid  of  the  previously  mentioned 
condition  of  the  tangency  in  the  plane  of  the  hodograph,  Karman 
seeks  an  infinite  acceleration  on  the  contour  of  the  NACA-4412 
profile,  which  has  no  infinite  cui^atures,  making  use  of  the  exper- 
imental distribution  of  the  velocity  on  the  contour.   In  the  plane 
of  the  hodograph,  Karman  obtains  the  point  of  tangency  of  the 
transform  of  the  contour  with  the  epicycloid  but  corresponding  to 
this  point  of  the  contour  of  the  profile,  infinite  acceleration 
does  not  occur  as  should  be  the  case.  Wo  satisfactory  explanation 
of  this  fact  contradicting  his  results  is  given. 
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This  noncorrespondence  is  explained  by  the  fact  that  the  tan- 
gency  in  the  plane  of  the  hodograph  can  occixr  simultaneously  with 
the  condition  S\|//S0  =  »  at  the  point  considered  on  the  flow  plane. 
An  example  of  this  is  the  Meyer  flow,  that  is,  the  flow  about  a 
certain  contour  for  which  the  entire  contour  coincides  in  the  plane 
of  the  hodograph  with  the  epicycloid  but  nevertheless,  at  the  points 
of  this  contour,  the  acceleration  is  finite  if  the  curvature  is 
finite. 

Moreover,  under  the  conditions  of  the  problem  of  the  flow  with 
local  supersonic  region,  the  tangency  in  the  plane  of  the  hodograph 
of  the  transform  of  the  contour  with  the  epicycloid  is  equivalent 
to  condition  (2.26)  at  the  point  of  the  contour  corresponding  to 
the  point  of  tangency,  because  at  this  point  of  the  contour  the 
acceleration  is  finite  if  the  curvature  of  the  contour  is  finite. 

In  recent  times,  investigations  have  appeared  on  exact  solu- 
tions of  the  problem  of  flow  with  local  supersonic  regions  where 
the  equations  of  S.  A.  Chaplygin,  in  the  plane  of  the  hodograph, 
are  used.  Notwithstanding  the  fact  that  in  these  papers  examples 
are  given  of  the  existence  of  potential  flows  with  local  supersonic 
region,  the  value  of  these  investigations  is  limited  by  the  fact 
that  with  change  in  the  Mach  number  of  the  oncoming  flow  there 
occurs  simultaneously  a  deformation  of  the  contour  of  the  body. 
Due  to  this  fact,  points  with  infinite  acceleration  appear  within 
the  supersonic  region,  which  on  the  basis  of  the  theorem  proved  in 
section  II,  part  2  indicates  an  infinite  curvature  of  the  contour 
at  these  points. 

From  the  foregoing  considerations,  it  is  evident  that  it  is  a 
mistake  to  associate  the  occurrence  of  such  a  limiting  line  of  the 
contour  with  the  impossibility  of  a  potential  flow  with  local 
supersonic  region  about  a  fixed  profile. 

The  investigation  of  the  principal  properties  of  a  flow  with 
local  supersonic  region  must,  in  the  opinion  of  the  author,  be 
carried  out  directly  in  the  flow  plane. 

In  the  work  described  herein,  all  results  were  obtained 
directly  from  investigations  in  the  flow  plane. 


Translated  by  S.  Eeiss 
National  Advisory  Committee 
for  Aeronautics 
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Figure  5. 
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Figure  8. 
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Figure  22. 


NACA    TM   No.     12 13 


33 
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